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University  of  Oklahoma ,  Norman,  Oklahoma 

Abstract 

A  theory  is  formulated  for  the  small -amplitude  free  vibration  of  thick, 
circular  cylindrical  shells  laminated  of  bimodulus  composite  materials,  which 
have  different  elastic  properties  depending  upon  whether  the  fiber-direction 
strain  is  tensile  or  compressive.  The  theory  used  is  the  dynamic,  shear 
deformable  (moderately  thick  shell)  analog  of  the  Sanders  best  first-approxi¬ 
mation  thin-shell  theory.  By  means  of  tracers,  the  analysis  can  be  reduced 
to  various  simpler  shell  theories,  namely  Love's  first  approximation,  and 
Donnell's  shallow-shell  theory.  As  an  example  of  the  application  of  the  theory, 
a  closed-form  solutionis  presented  for  a  freely  supported  panel  or  complete 
shell.  To  validate  the  analysis,  numerical  results  are  compared  with  existing 
results  for  various  special  cases.  Also,  the  effects  of  the  various  shell 
theories,  thickness  shear  flexibility,  and  bimodulus  action  are  investigated. 
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Nomenclature 


A..  ■  stretching  stiffnesses 

*  J 

B--  3  bending-stretching  coupling  stiffnesses 

*  J 

b  3  circumferential  arc  width  of  panel 
[C]  3  matrix  of  coefficients 

C^.C^  3  shell -theory  tracers 

~Z1  s  VR 

D,._.  3  bending  stiffnesses 

d  3  32(  )/3X3y 

h  3  total  shell  thickness 

I  3  rotatory  inertia  coefficient  (per  unit  midsurface  area) 

K4,K5  3  thickness-shear  correction  factors 
[L]  3  matrix  of  linear  differential  operators 

3  stress  couples  and  in-surface  stress  resultants 
m,n  3  axial  and  circumferential  mode  numbers 
P  *  normal  inertia  coefficient  (per  unit  midsurface  area) 

Qi  3  thickness-shear  stress  resultants 

Qi jkji  "  Plane~stress  reduced  stiffness  coefficients 

R  3  radius  of  shell  midsurface 

Sjj  3  thickness-shear  stiffnesses 

t  3  time 

u,v,w  3  axial,  circumferential ,  and  radial  midsurface  displacements 
U,V,W  3  amplitudes  of  u.v.w 

x,y,z  3  position  coordinates  In  axial,  circumferential,  and  outward  normal 
directions 

X,Y  3  amplitudes  of  and  \i>2 
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znx,zny 

Zx’Zy 

а,  6 

б,  A 

'l»l  3^2 


P 

ai 

ai 

(  ). 


xy 


neutral -surface  positions  forex*0  and  ey*0*  respectively 

■  V’  V 

=  niir/L,  nir/b  (for  panel)  or  n/R  (for  complete  cylinder) 

*  generalized  displacement  vectors  defined  in  Eqs.  (8)  and  (15) 

*  bending  slopes  in  xz  and  yz  planes 

*  strain  components  at  arbitrary  location  (x,y,z) 

*  midsurface  strain  components 

*  curvature-change  components 
=  material  density 

■  stress  components 

*  natural  frequency 

*  a2(  )/axay 
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Introduction 

Analysis  of  composite-material  shell-type  structures  has  been  of  con¬ 
siderable  research  interest  in  the  past  two  decades.  Earlier  work  has  been 
summarized  in  a  monograph  by  Leissa1  and  in  two  more  recent  survey  papers  ’  . 

Due  to  its  importance,  the  finite-length  circular  cylindrical  shell  con¬ 
figuration  has  received  the  most  attention.  Thin-shell  analyses  of  composite- 
material  shells  of  this  configuration  have  been  performed  by  Dong4,  Abhat  and 
Wilcox5  and  Fortier  and  Rossettos6  using  Donnell’s  shallow-shell  theory7,  by 
Bert,  Baker,  and  Egle8,  Stavsky  and  loewy9,  Shivakumar  and  Krishna  Murty10, 
and  Greenberg  and  Stavsky11  using  Love's  first  approximation  theory12,  and  by 
Padovan13  using  Novozhllov's  higher-order  theory14.  In  the  present  work, 
tracer  coefficients  are  introduced  to  consider  Sanders'  best  first-approximation 
theory15.  Love's  first  approximation12,  Loo's16,  Morley's17,  and  Donnell's7 
theories. 
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The  thickness  shear  moduli  of  composite  materials  are  typically  an  order 
of  magnitude  lower  than  their  in-plane  elastic  moduli.  This  is  in  contrast 
to  isotropic  materials,  which  have  shear  to  elastic  moduli  ratios  of  approxi¬ 
mately  40%.  Thus,  composite-material  shells  would  be  expected  to  be  affected 
by  thickness-shear  flexibility  to  a  much  greater  extent  than  isotropic  shells 
having  the  same  geometric  parameters.  Previous  research  on  vibration  of 
laminated,  shear  deformable  shells  includes  that  of  Refs.  18  and  19,  using  the 
shear  deformable  analog  of  Donnell's  shallow  shell  theory,  and  Refs.  20  and  21, 
using  the  shear  deformable  version  of  Love's  first  approximation.  In  some 
recent  work22,  the  shear  deformable  analog  of  Sanders'  theory  was  developed 
and  it  was  shown,  using  tracers,  that  this  theory  can  be  reduced  to  two  simpler 
shear  deformable  theories  as  special  cases:  (1)  Love's  first  approximation 
and  Loo's,  and  (2)  Morley's  and  Donnell's.  This  approach  is  extended  to  shell 
vibration  in  the  present  paper, 

Certain  fiber- reinforced  materials,  especially  those  with  yery  soft 
matrices  (e.g.,  cord-rubber),  exhibit  quite  different  elastic  behavior  depend- 

23-25 

Ing  upon  whether  the  fiber-direction  strain  is  tensile  or  compressive 
As  a  first  approximation,  the  stress-strain  behavior  of  such  materials  is 
usually  modeled  as  being  bilinear,  with  different  elastic  coefficients 
depending  upon  the  sign  of  the  fiber-direction  strain.  Thus,  these  materials 
are  called  blmodulus  composite  materials.  The  fiber-governed,  symmetric- 
compliance  material  model  proposed  In  Ref.  26  has  been  shown  to  agree  well 
with  experimental  data  obtained  under  in-plane  loading  of  several  materials 
with  drastically  different  properties  in  tension  and  compression.  To  the 
best  of  the  present  Investigators'  knowledge,  all  of  the  previous  investigations 
of  laminated  bimodulus  shells,  as  reviewed  in  Refs.  22  and  27,  have  Involved 
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the  static  behavior  of  thin  cylindrical  shells  under  mechanical  or  thermal 
loading.  However,  free  vibration  of  cross*ply  laminated  plates  was  considered 
In  Ref.  28. 

Formulation 

Laminate  Action 

To  explain  the  effect  of  bimodulus  action  on  the  shell  behavior,  we 
discuss  first  a  flat  plate.  If  a  single-layer  plate  or  a  laminate  symmetrically 
laminated  about  its  midplane  are  made  of  ordinary  (not  bimodulus)  materials, 
there  is  no  coupling  between  bending  and  stretching.  However,  a  laminate 
which  is  not  symmetrically  laminated  does  exhibit  bending-stretching  coupling. 

Now  consider  a  flat  plate  consisting  of  either  a  single  layer  of  bimodulus 
material  or  symmetrically  laminated  of  bimodulus  materials.  When  the  plate 
is  subjected  to  bending,  the  different  elastic  moduli  in  tension  and  compres¬ 
sion  cause  a  shift  in  the  neutral  surface  away  from  the  geometric  midplane 
toward  the  tension  side  of  the  bend  if  the  tensile  properties  exceed  the  com¬ 
pressive  ones.  Consequently,  this  plate  displays  bending-stretching  coupling, 
analogous  to  a  plate  unsymmetrically  laminated  of  ordinary  material  as  des¬ 
cribed  in  the  preceedlng  paragraph.  Therefore,  in  a  laminated  bimodulus- 
material  shell,  bending-stretching  coupling  can  be  induced  by  the  combination 
of  bimodulus-material  action,  laminate  configuration,  and  shell  geometry 
(curvature). 


One  can  write  the  laminate  constitutive  relations  as 


(Q^  »  CS1j]  {ej}  ( 1  »j*4,5) 
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Here,  the  respective  In-plane,  bending-stretching,  bending  or  twisting,  and 
thickness-shear  stiffnesses  are  defined  as 

h/2 

<WV  ■  j  dz  a.J-1.2,6) 

-h/2  {2) 

h/2 

stj  ■  \  KiKj  w dz  (1-J’4i5> 

-h/2 

Here,  denotes  the  plane-stress  reduced  stiffness  defined  by 

{ai}  *  CQijkt:i{ej}  (3) 

and  subscripts  1j  refer  to  the  position  In  the  stress-strain-relation 
array,  k  refers  to  the  sign  of  the  fiber-direction  strain  (1  ^  +,  2  ^  -), 
and  i  Is  the  layer  number. 

Shear  Deformable  Shell  Theory 

The  coordinate  system  used  here  Is  shown  in  Fig,  1.  Then  the  equations 
of  motion  In  the  absence  of  body  forces  and  body  moments  can  be  written  as 


Figure  1.  Shell  geometry 


7 


I 


v 


*  "e.y  '  <V2R>Vy  *tt 

"6>x  +  N2,y  *  'V™.  +  «V2R>M6,x  -  *  <^R>*2,tt 

Qs.x  *  \,y  -  (VR)  '  Pw-tt 
"i,x  *  Ms,y  ‘  «,  '  **!.«  * 

M6,x  *  M2,y  '  \  ■  '*2,tt  * 


(4) 


Here,  it  has  been  assumed  that  the  mass  distribution  is  symmetric  about  the 
middle  surface.  The  normal  and  rotatory  inertias  are 


(P.I)  - 


h/2 

-h/2 


(l,z2)p  dz 


(5) 


The  first  terms  on  the  right  side  of  each  of  Eqs.  (4)  are  the  same  as 
for  plates  in  Mindlin's  shear  deformable  theory29.  The  second  terms  on  the 
right  side  of  all  except  the  third  equation  of  set  (4)  are  the  inertia-coupling 
terms  due  to  shell  curvature,  first  introduced  by  Mirsky  and  Herrmann30. 

The  kinematics  of  deformation  can  be  expressed  as 


e?  *  u. 


1  *  v,„  +  (w/R) 

’y 


1  “’x  ’  "2  y 

e4  »  «J>2  +  w,  -  (C1/R)v  ;  eg  ■  +  w. 


;  e°  s  U.  +  v, 

6  *y  x 


(6) 


Ki  *  *i,x 


^.y  ’  k6  3  ^l.y  +  *2>x  +  (C2/2R)(v,x*u,  ) 


The  coefficients  appearing  In  Eqs.  (4)  and  (6)  are  shell-theory 
tracers  which  take  on  values  listed  In  Table  1.  Note  that  in  contrast  to 
the  case  of  thln-shell  theory.  Loo's  approximation  offers  no  simplification 
over  Love's  first-approximation  theory  and  Donnell's  theory  coincides  with 
Morley's. 
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Table  1.  Shell -theory  tracers  and  their  values 


Theory  (thin- she 11  theory 
generalized  to  shear  deformable  theory) 

Ci 

c2 

Sanders' 15 

1 

1 

Love's  first  approximation12  and  Loo's16 

1 

0 

Morley's17  and  Donnell's7 

0 

0 

Substituting  Eqs.  (1)  and  (6)  into  Eqs.  (4),  one  obtains  the  following 


operator  equation 

[L]{<5}  =  0 

(7) 

where 

{6}  =  {u.v.w,^,^}1 

(8) 

and  [L]  is  a  symmetric  matrix  of  linear  differential  operators  listed  in 
Appendix  A. 

Criteria  for  Homogeneity  Along  the  Middle  Surface 

In  deriving  Eq.  (7),  we  tacitly  assumed  that  the  laminate  stiffnesses 
(A.  .,B.  .,D..,S. .)  are  all  independent  of  coordinates  (x,y)  on  the  middle 

I  J  I  J  I  J  •  J 

surface.  However,  in  view  of  the  bimodulus  nature  of  the  materials  comprising 
the  laminate,  these  stiffnesses  depend  upon  the  fiber-direction  neutral- 
surface  positions  associated  with  the  respective  layers  (i.e.,  znx  for  a 
single  layer  with  axially  oriented  fibers,  and  znx  and  zny  for  a  cross-ply 
laminate). 

Thus,  for  layers  having  the  fibers  oriented  axially,  the  associated 
fiber-direction  neutral -surface  position  is  determined  by 

ei  8  e?  +  znxKi  =  0 
or 

znx  “  "  CVK\  8  '  u,x/4,i,x  8  constant  (9) 
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Similarly,  for  layers  having  the  fibers  oriented  circumferentially 

£_  =  £3  +  Z  K.  *  0 
2  2  ny  2 

or 

znv  =  “  eVs  =  -(v»w  +  =  constant  (10) 

Closed-Form  Solution 

We  seek  a  solution  which  satisfies  the  governing  operator  Eq.  (7),  the 
homogeneity  relations  (9)  and  (10),  and  appropriate  boundary  conditions. 

A  closed-form  solution  is  presented  for  a  shell  that  is  freely  sup¬ 
ported  (simply  supported  without  in-surface  restraint;  SS3  in  Hoff's 
notation31)  along  its  curved  edges: 

Ni(0,y,t)  =  Ni(L,y,t)  =  M^O.y.t)  =  Mi(L,y,t)  =  0 

w(0,y,t)  =  w(L,y,t)  *  v(0,y,t)  =  v(L,y,t)  =  0  (11) 

02(O,y,t)  *  iP2(L,y,t)  =  o 

The  shell  may  be  either  circumferentially  complete  (closed),  which 
requires  that  the  circumferential  mode  shapes  be  periodic  in  position  (i.e., 
trigonometric  functions),  or  a  cylindrically  curved  panel  (open  shell)  freely 
supported  along  the  generators,  which  requires  that  the  following  additional 
boundary  conditions  be  satisfied: 

Mx.O.t)  *  Mx.b.t)  =  M2(x,0,t)  =  M2(x,b,t)  =  0 

w(x,0,t)  *  w(x,b,t)  =  u(x,0,t)  *  u(x,b,t)  *  0  (12) 

i^j(x,0,t)  »  ^(x.b.t)  =  o 

Under  these  conditions,  the  solution  to  Eq.  (7)  is  of  the  form 


i 
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u  =  U  cos  ax  sin  By  cos  wt 

v  =  V  sin  ax  cos  By  cos  u>t 

w  *  W  sin  ax  sin  By  cos  ait  (13) 

*  X  cos  ax  sin  By  cos  ait 

i|>2  *  Y  sin  ax  cos  By  cos  wt 

where  a=mir/L  and  B=nir/b  for  a  panel  and  n/R  for  a  complete  cylinder. 

Substitution  of  Eqs.  (13)  into  Eq.  (7)  leads  to  the  following  homo¬ 
geneous  algebraic  system 

[C]{A>  =  {0}  (14) 

where  [C]  is  a  matrix  with  elements  depending  upon  A..,  B..,  D..,  S..,  P,  I, 

IJ  1 J  ^  J 

a,  and  s,  and 

{A}  =  {U,V,W,X,Y}T  (15) 


Numerical  Results 

Before  proceeding  to  the  calculations  for  shells  laminated  of  bimodulus 
materials,  we  present  some  results  for  both  circumferentially  complete 
cylinders  and  cylindrically  curved  panels  laminated  of  ordinary  composite 
materials. 

The  first  example  is  for  a  thin-walled  complete  cylinder  cross-plied 
of  boron-epoxy.  The  geometric  parameters  and  material  properties  used  are 
listed  in  Table  2.  The  agreement  between  the  natural  frequencies  for  numer¬ 
ous  modes  computed  by  use  of  the  thin-shell  theory  of  Appendix  B  (special¬ 
ized  to  Love's  first  approximation)  and  the  Love’s  first-approximation 
analysis  of  Ref.  8  is  excellent,  as  shown  by  comparison  of  columns  3  and  4 
in  Table  3.  The  inclusion  of  thickness  shear  deformation  sometimes 


1 


n 


decreases  the  natural  frequencies,  but  in  other  instances  increases  them, 
as  can  be  seen  in  column  5  of  Table  3.  For  this  thin  shell,  the  effect 
of  rotatory  inertia  is  almost  negligible  (column  6). 

The  second  example  is  for  a  cyl indrically  curved  panel  with  thickness 
shear  flexibility  included,  using  numerical  input  data  from  Refs.  18  and 
19  as  listed  in  Table  4.  In  order  to  make  a  fair  comparison  with  the 
results  of  Refs.  18  and  19,  it  is  necessary  to  use  the  shear  deformable 

version  of  Donnell 's  theory  and  to  neglect  rotatory  inertia  and  the  tangen- 

32 

tial  inertias  (Pu,tt  and  Pv,tt)  as  suggested  by  Vlasov  .  Dimensionless 
natural  frequencies  for  various  modes,  with  and  without  bending-stretching 
coupling  (B.,.)  suppressed,  are  listed  in  Tables  5  and  6.  The  lack  of  good 

*  J 

agreement  in  Table  5  is  difficult  to  explain,  especially  in  view  of  the 
good  agreement  in  Table  6. 

In  studying  the  above  results,  one  can  discern  these  general  trends: 

1.  For  a  complete  shell  at  a  fixed  axial  wave  number,  the  natural 

frequencies  first  decrease  then  increase  with  increasing  circumferential 

33 

wave  number  (n).  As  was  first  explained  by  Arnold  and  Warburton  ,  this 
is  due  to  the  decrease  in  membrane  strain  energy  and  the  increase  in  bending 
strain  energy  as  n  is  increased. 

2.  The  inclusion  of  thickness  shear  deformation  and  rotatory  inertia 

does  not  always  decrease  the  frequency.  A  similar  observation  may  be  seen 

34 

in  the  recent  work  of  Greenberg  and  Stavsky  . 

3.  Bending-stretching  coupling  (B^.^0)  always  lowers  the  natural 
frequency  compared  to  the  value  for  the  uncoupled  case  (B^.  =  0). 

4.  If  the  frequencies  are  normalized  with  the  panel  width  (b),  as  in 
Table  6,  increasing  the  aspect  ratio  (L/b)  lowers  the  frequency  values, 
just  as  it  does  in  the  case  of  a  flat  panel  (plate). 
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To  illustrate  the  effects  of  bimodulus  composite-material  action, 

aramid-rubber  properties. determined  primarily  from  the  experiments  of 
24 

Patel  et  al.  by  use  of  the  material  model  of  Ref.  26,  are  used.  Some 
details  of  the  estimation  of  certain  properties  not  directly  measured  in 
Ref.  24  are  given  in  Ref.  28.  The  complete  set  of  properties  used  is 
listed  in  Table  7.  This  material  has  the  most  drastic  difference  between 
tensile  and  compressive  properties  known  to  the  present  investigators; 
thus,  it  is  definitely  a  worst  case. 

Although  the  present  problem  is  a  linear  one  in  each  portion  of  a 
cycle  of  vibration,  it  actually  is  different  in  the  two  portions  of  a  cycle. 
This  phenomenon  was  explained  for  the  vibration  of  single-layer  and  two- 
layer,  cross-ply  plates  in  Ref.  28.  Here,  for  completeness, we  discuss  the 
cross-ply  case:  side  cross-sectional  views  are  shown  in  Fig.  2.  The  inner 
layer  is  oriented  circumferentially  (y  direction)  and  the  outer  layer  is 
oriented  axially  (x  direction).  During  the  first  portion  of  a  cycle,  con¬ 
sider  the  shell  to  be  deflected  as  shown  in  Fig.  2(a).  Then  the  neutral 
surface  for  ex  falls  outside  of  the  interface  within  the  axially  oriented 
layer,  while  the  neutral  surface  for  Cy  lies  inside  the  interface,  i.e., 
within  the  circumferential  layer.  In  the  latter  portion  of  the  vibration 
cycle.  Fig.  2(b),  the  ex  neutral  surface  falls  outside  the  axial  layer,  and 
the  Cy  neutral  surface  lies  outside  of  the  circumferential  layer.  Therefore, 
in  the  latter  portion  of  the  cycle,  compressive  properties  are  used  for 
the  entire  axial  layer  and  tensile  properties  for  the  entire  circumferential 
layer.  From  the  foregoing  considerations  for  a  two-layer,  cross-ply  shell, 
it  is  clear  that  the  stiffnesses  acting  in  the  two  portions  of  a  cycle, 

(a)  and  (b)  in  Fig.  2,  are  different  and  thus  the  associated  frequencies 


(denoted  by  on  and  u>2,  respectively)  are  also  different.  The  corresponding 
time  intervals  over  which  the  two  portions  take  place  are  ir/o>i  and  it /w2. 
Thus,  the  total  period  for  a  complete  cycle  is  (u/ui)  +  (ir/ui2).  The  average 
frequency  (u>)  over  the  entire  cycle  is  2ir  divided  by  the  total  period;  then 

w  1  *  (l/2)(u>i  +  <i>2  )  (16) 

As  discussed  above,  the  neutral-surface  positions  are  different  during 
the  two  portions  of  a  cycle;  numerical  results  are  listed  in  Table  8.  The 
average  frequency,  as  calculated  from  Eq.  (16),  is  plotted  in  Fig.  3  versus 
aspect  ratio  for  various  values  of  R/h.  The  dashed  line  shows  the  analogous 
results  from  Ref.  28.  for  flat  plates. 

As  was  mentioned  in  Ref.  28,  it  can  be  shown  that  in  the  bilinear 
analysis  presented,  energy  is  conserved  at  the  junction  between  the  two 
portions  of  the  vibrational  cycle. 
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Concluding  Remarks 

Two  analyses  have  been  presented  for  prediction  of  the  free-vibrational 
behavior  of  freely  supported,  cylindrical ly  curved  panels  or  complete 
cylindrical  shells  cross-ply  laminated  of  bimodulus  composite  materials. 

One  analysis  is  for  thin  shells  and  the  other  is  for  moderately  thick 
shells  with  thickness  shear  deformation  and  rotatory  inertia  included. 

When  reduced  to  the  special  case  of  ordinary  material,  the  solutions  yielded 
results  which  agreed  well  with  various  existing  numerical  results.  The 
peculiarities  of  bimodulus  shells  have  been  discussed  and  the  effects  of  a 
limited  number  of  parameters  investigated.  These  results  may  be  used  to 
validate  numerical  analyses,  such  as  those  based  on  the  finite-element 
method . 
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Appendix  A:  linear  Differential  Operators  Appearing  In  Eg.  (7) 
1-11  *  Aud2  +  (^66  "  ^2®66  +  J  ^2^66^y  '  Pdt 

L12  *  (A12  +  ^66  "  J 

L 1 3  *  (A12/R)dx  ;  Li4  *  Bud2  +  (®66  "  \  C2D66)d2  " 

L i 5  *  (B12  +  B66  -  C2066)dxdy 

L22  *  (*66  +  C2B66  +  ^  C2D66)dx  +  A22d2  -  C2S44  -  Pd^ 

L23  *  (R  lft22  +  C 1S44 )dy  ;  L24  *  (Bu  +  B66  +  j  C2066)dxd^ 

L25  *  (066  +  J  ^2^66^x  +  ®22^y  +  ^1^44  “  (I/R)d^ 

L33  *  ■  Sssd2  “  S44d^  +  (A22/R2)  +  Pd2 


(R  lBi2  -  S5 5 )dx  ;  L35  * 

(R  -  S44 )dy 

Dudx  +  066<ly  -  S55  -  Id2 

;  L45  *  (D12  +  D66)dxdy 

D66dx  +  ^22dy  ”  $44  -  Id^ 

;  »  C4/R  ;  dx  -  »(  )/3x 

Appendix  B: 

Th In-Shell  Theory 

The  theory  developed  here  Is  an  extension  of  the  static  bimodulus- 
composite  shell  analysis  presented  In  Ref.  27.  The  shell  theory  used  is  a 
generalized  first-approximation  one  that  can  be  reduced  to  the  theories  of 
Sanders,  Love's  first-approximation.  Loo,  Morley,  and  Donnell  by  use  of  the 
tracers  listed  In  Table  B-l.  It  is  emphasized  that  these  are  different  tracers 
than  those  listed  In  Table  1. 

The  shell  constitutive  equations  are  the  same  as  given  In  Eqs.  (l)-(3), 
except  now  the  thickness  shear  quantities  (Q4  and  Qg)  are  not  needed.  The 
appropriate  equations  of  motion  are 


i 


"i.x  *  \.y  •  <c;/R)M6.y  ■  Pu-tt 


>6.x  *  "2.y  *  <c;/r!m6.,  +  <c;/R)«a,y  '  p»-tt 
'«.«x  *  “s.xy  *  M2,yy  '  <VR>  *  ”»•« 


The  kinematic  relations  for  the  middle-surface  strains  are  the  same  as 
in  Eqs.  (6).  Since  now  the  thickness  shear  strains  (e4  and  e5)  are 
identically  zero,  the  curvatures  are  given  by 


‘1  *  •  "'xx  1  "2  ■  •  w’yy  *  <C;'R>v'y 

*6 =  -  ^-xy  *  <q'Rw.x  •  (c;/R>u-y 


Substituting  the  expressions  for  the  middle-surface  strains  and 
curvatures  into  the  constitutive  relations  and  the  latter  into  Eqs.  (Bl), 
one  obtains  the  same  form  of  operator  expression  as  Eq.  (7)  except  that 
now  the  elements  of  the  {5}  vector  are  only  u,  v,  w.  The  linear  differ¬ 
ential  operators  are  now 

L 1 1  *  And*  +  [A66  -  2C;b66R"1  +  (c;)2D66R’2]dJ  -  Pd2 
L 12  -  [A  12  +  A66  +  C2B12R" 1  +  (C!)  -  CUB66R_1  -  C^C;D66R'2]dxdy 
L 1 3  *  -  &ndx  -  (B 12  +  2B66  -  2C4D66R  )^xd2  +  Al2R  *dx  *  L31 
L2i  -  [A  12  +  A66  +  CiB12R_1  +  (C^  -  C;)B66R_1  -  C^C;D66R'2]dxdy 
l22  *  CA66  +  2C’3B66R‘1  +  (C^)2D66R-2]dJ 

+  [A22  +  (C'l  +  Ci)B22R“1  +  CiC£D22R’2]dy  -  P<i\ 

L23  *  "  CB 12  +  2Bgg  +  C1DJ2A  +  C^DggR  l]d2dy 

-  (B22  +  C'1022R"1)dy  +  (A22  +  CiB22R"1)R’ldy 
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1-32  *  “  CB 12  +  2B6g  +  C2Dl2R  +  2C3DggR 

-  (B22  +  C£D22R-1)dy  +  (A22  +  C^R'V1^ 

L 33  *  +  ^°12  +  2D66)d2d2  +  D22dy 

-  2B12R'ld2  -  2B22R”1dy  +  A22R’2  +  Pd2  (B3) 

Note  that  the  operators  Lr$  are  symmetric  for  all  of  the  theories  except 
Morley's. 

The  criteria  for  the  neutral -surface  locations  (znx  and  zny)  to  be 
constant  are  now 


znx  =  u’x/w*xx  *  const‘ 

Zny  *  ^V’y  +  '  C2R"ly*y)  *  const. 


(B4) 


The  appropriate  boundary  conditions  for  closed  cylindrical  shells  and 
cylindrical ly  curved  panels  are  the  same  as  Eqs.  (11)  or  Eqs.  (12)  except 
that  here  the  conditions  on  and  ii2  are  not  needed  (since  ifi  and  <f>2  do 
not  appear  In  this  thln-shell  theory).  Similarly,  only  the  first  three 
expressions  in  Eqs.  (13)  and  the  first  three  elements  In  Eqs.  (15)  are 
required. 
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Table  2.  Geometric  Parameters  and  Material -Property 
Data  for  Example  I  (Ref.  8) 


Geometric  Parameters 
Mall  thickness  h.  In  (mm) 

Middle-surface  radius  R,  In  (mm) 

Length  Lt  In  (mm) 

Material  Properties  (Boron-Epoxy) 
Fiber-direction  Young's  modulus,  psl  (GPa) 
Transverse  Young's  modulus,  psl  (GPa) 

Major  Poisson's  ratio,  dimensionless 
In-plane  and  thickness  shear  moduli,  psl  (GPa)*2 
Specific  gravity,  dimensionless 


0.02  (0.51) 

2.481  (63.0) 

31.5  (800) 

31.0  x  106  (214) 
2.7  x  106  (18.6) 
0.28 

0.75  x  106  (5.17) 
2.05 


a  The  minor  Poisson's  ratio  Is  assumed  to  be  given  by  the  reciprocal 
relation. 


Table  3.  Natural  Frequencies  Associated  with  Various  Vibrational  Modes  for 
a  Complete  Cylindrical  Shell  of  Boron-Epoxy  (Example  I) 


No.  of  axial 
half  waves. 

No.  of  circum. 
full  waves, 

Natural  Frequencies, 

Hz 

Ref.  8a 
TST 

Present 

m 

n 

TST 

SOT 

SOT  wRI 

1 

1 

532 

534 

533 

534 

2 

235 

241 

240 

246 

3 

253 

260 

256 

257 

4 

444 

450 

418 

412 

5 

714 

721 

751 

823 

6 

1047 

1056 

1057 

1113 

7 

1442 

1453 

1460 

1536 

8 

1897 

1911 

1941 

2028 

2 

1 

1287 

1290 

1290 

1290 

2 

676 

683 

682 

682 

3 

443 

454 

445 

444 

4 

497 

509 

480 

505 

5 

728 

739 

753 

830 

6 

1051 

1063 

1083 

1135 

7 

1442 

1456 

1489 

1511 

8 

1897 

1913 

1932 

2020 

a  The  numerical  values  listed  here  were  taken  directly  from  the 
computer  printout,  rather  than  from  the  less  accurate  graphical  pre 
sentation  of  Ref.  8,  which  was  not  plotted  accurately. 


TST:  thin-shell  theory 
SOT:  shear  deformable  theory 

SOT  wRI:  shear  deformable  theory  with  rotatory  inertia 


l 
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Table  4.  Geometric  Parameters  and  Material -Property 
Data  for  Example  II 


Geometric  Parameters 

Ex.  IIA 
(Ref.  18) 

Ex.  I IB 
(Ref.  19) 

Panel  aspect  ratio  (L/b) 

1 

1  to  5 

Middle-surface  radius/thickness  (R/h) 

31.25  to  312.5 

40 

Material  Properties 

Fiber-direction  Young’s  modulus 

40  Et 

25  Et 

Transverse  Young's  modulus 

Et 

et 

Major  Poisson's  ratio*2 

0.25 

0.25 

In-plane  and  longitudinal-thickness 
shear  moduli 

et 

0.5  Et 

Transverse-thickness  shear  modulus 

0.5  Et 

0.2  Et 

a  The  minor  Poisson's  ratio  is  assumed  to  be  given  by  the 

reciprocal  relation. 
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Table  5.  Dimensionless  Natural  Frequencies  uib^P/Eyh3)15  Associated  with  Various 
Vibrational  Modes  for  Cylindrically  Curved  Panels  (Example  IIA) 


No.  of  axial 
half  waves 
m 

No.  of  circ. 
half  waves 
n 

1 

R/h 

w/o  Coupling 

Ref.  18 

Present 

w/Coupling 

Ref.  18 

<v°> 

Present 

1 

1 

50 

312.5 

24.76 

20.59 

19.23 

14.08 

2 

1 

50 

312.5 

53.44 

55.40 

- 

- 

1 

2 

50 

312.5 

60.43 

53.58 

- 

- 

2 

2 

50 

312.5 

75.86 

75.51 

- 

- 

1 

1 

50 

156.25 

36.00 

24.48 

32.86 

19.32 

2 

1 

50 

156.25 

54.91 

- 

- 

1 

2 

50 

156.25 

78.74 

53.95 

- 

- 

2 

2 

50 

156.25 

80.34 

76.66 

- 

- 

1 

1 

10 

62.5 

22.00 

16.91 

18.57 

11.34 

2 

1 

10 

62.5 

35.06 

38.30 

- 

- 

1 

2 

10 

62.5 

45.01 

38.19 

- 

- 

2 

2 

10 

62.5 

48.46 

51.70 

- 

- 

1 

1 

10 

31.25 

34.86 

17.11 

32.38 

11.64 

2 

1 

10 

31.25 

37.27 

38.63 

- 

- 

1 

2 

10 

31.25 

67.63 

38.21 

- 

- 

2 

2 

10 

31.25 

55.21 

51.77 

- 

- 
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Table  6.  Dimensionless  Fundamental  Frequencies  ojb2(P/Eyh3)ls 
for  Cylindrical ly  Curved  Panels  (Example  I IB ) 


Panel  Aspect 
Ratio  L/b 

Without  Coupling 

With  Coupling 

Ref.  19 

Present 

Ref.  19 

Present 

1 

16.11 

15.89 

11.71 

11.65 

2 

10.94 

10.79 

7.35 

7.37 

3 

10.27 

10.10 

6.58 

6.59 

4 

10.07 

9.91 

6.32 

6.33 

5 

9.99 

9.83 

6.22 

6.21 

Table  7.  Material  Properties  for  Aramid-Rubber 
Composite  Material* 

Bimodulus 

Property 

k  *  1 

II 

Fiber-direction  Young's  modulus,  GPa 

3.58 

0.0120 

Transverse  Young's  modulus,  GPa 

0.00909 

0.0120 

Major  Poisson's  ratio,  dimensionless 

0.416 

0.205 

Minor  Poisson's  ratio,  dimensionless 

0.00106 

0.205 

In-plane  and  longitudinal-thickness 
shear  moduli,  GPa 

0.00370 

0.00370 

Transverse-thickness  shear  modulus,  GPa 

0.00290 

0.00499 

Specific  gravity,  dimensionless 

0 

.970 

a  Fiber-direction  tension  is  denoted  by  k  =  1,  and  fiber- 
direction  compression  by  k  =  2. 
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Table  8.  Dimensionless  Neutral -Surface  Locations  in  the  First  and 

Second  Portions  of  a  Cycle  for  Two-Layer,  Cross-Ply  Cylindrical 
Panels  Constructed  of  Aramid-Rubber  (Sanders-Type  Theory  with 
Thickness  Shear  Deformation;  b/h  =  10) 


Aspect 

ratio 

a/b 

First  Portion 

Second 

Portion 

*.n> 

z  (’> 

_ y _ 

z  <2> 

X 

2,<2) 

R/h  =  10  (Highly  Curved) 

0.5 

0.4405 

-0.0987 

-0.0272 

0.3824 

0.7 

0.4343 

-0.0736 

-0.0374 

0.4108 

1.0 

0.4234 

-0.0521 

-0.0546 

0.4283 

1.4 

0.4069 

-0.0368 

-0.0768 

0.4378 

2.0 

0.3796 

-0.0250 

-0.1101 

0.4433 

R/h  =  20  (Moderately  Curved) 

0.5 

0.4430 

-0.0327 

-0.0226 

0.4033 

0.7 

0.4387 

-0.0618 

-0.0316 

0.4220 

1.0 

0.4312 

-0.0438 

-0.0459 

0.4338 

1.4 

0.4195 

-0.0311 

-0.0648 

0.4401 

2.0 

0.3995 

-0.0213 

-0.0930 

0.4439 

R/h  =  50  (Slightly  Curved) 

0.5 

0.4446 

-0.0719 

-0.0193 

0.4161 

0.7 

0.4418 

-0.0539 

-0.0272 

0.4289 

1.0 

0.4361 

-0.0384 

-0.0394 

0.4371 

1.4 

0.4276 

-0.0273 

-0.0558 

0.4414 
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Table  B-l  Thin-shell -theory 

tracers 

and 

their 

values 

Theory 

Cl 

Cl 

C3 

c; 

Sanders 

1 

1 

3/2 

1/2 

Love’s  first  approximation 

1 

1 

1 

0 

Loo 

1 

0 

1 

0 

Morley 

1 

0 

0 

0 

Donnell 

0 

0 

0 

0 

8 


AXIS  OF  SHELL- 


v° 

=  0 


(a)  First  portion  of  cycle 


(b)  Second  portion  of  cycle 


Bimodulus  action  during  the  two  portions  of  motion  of  a  two-layer, 
cross-ply  shell  in  the  fundamental  mode  of  vibration.  Inside  layer 
is  in  y  direction  (90°);  outside  layer  is  in  x  direction  (0°).  Shaded 
portions  are  in  tension  in  the  respective  fiber  directions. 


ASPECT  RATIO,  L/b 


Fig.  3  Variation  of  fundamental  frequency  with  aspect  ratio  for 
two-layer,  cross-ply  curved  panels  of  aramid-rubber. 
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